Local Projections of Low-Momentum Potentials 
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Nuclear interactions evolved via renormalization group methods to lower resolution become in- 
creasingly non-local (off-diagonal in coordinate space) as they are softened. This inhibits both the 
development of intuition about the interactions and their use with some methods for solving the 
quantum many-body problem. By applying "local projections", a softened interaction can be re- 
duced to a local effective interaction plus a non-local residual interaction. At the two-body level, 
a local projection after similarity renormalization group (SRG) evolution manifests the elimination 
of short-range repulsive cores and the flow toward universal low- momentum interactions. The SRG 
residual interaction is found to be relatively weak at low energy, which motivates a perturbative 
treatment. 

PACS numbers: 21.30.-x,05.10.Cc,13.75.Cs 



I. INTRODUCTION 

Nuclear two-body interactions found in textbooks are 
readily visualized because they are usually local; that 
is, they are radially diagonal in coordinate representa- 
tion (e.g., see Refs. [1-3]). Plots of the central part of 
the potential, for example, exhibit long-range attraction 
from one-pion exchange, moderate mid-range attraction, 
and strong short-range repulsion. The long-range behav- 
ior is established from general considerations of sponta- 
neously broken chiral symmetry [4] but the short-range 
features are tied to the constraint of locality. That is, 
a repulsive core is the inevitable consequence of fitting a 
local nucleon-nucleon (NN) interaction to S-wave scatter- 
ing phase shifts beyond energies where they turn repul- 
sive. While short-range locality is not a physical require- 
ment [5], the associated short-range correlations (SRC) 
induced by the repulsive core provide intuition for nu- 
clear observables, such as those associated with high- 
momentum-transfer probes [6]. 

Local potentials are an advantage or even a neces- 
sity for some methods of solving the nuclear many-body 
problem. In particular, current implementations of quan- 
tum Monte Carlo, such as Green's function Monte Carlo 
(GFMC) [7, 8] and auxiliary field diffusion Monte Carlo 
(AFDMC) [9] require interactions that are expressed in 
an operator basis multiplying local functions of internu- 
cleon distances (e.g., Argonne vi^ [10])- On the other 
hand, the repulsive core makes such potentials ill-suited 
for basis expansions used in configuration interaction 
(CI) and coupled cluster (CC) methods because the ba- 
sis sizes required for convergence grow too large. In- 
stead one turns to softened nuclear potentials, which sup- 
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press coupling between low- and high-momentum com- 
ponents and thereby exhibit greatly improved conver- 
gence. Examples of soft NN potentials include chiral ef- 
fective field theory potentials with cutoffs of order 500 
to 600 MeV [4, 11, 12] and renormalization group (RG) 
evolved interactions [13]. 

These low-momentum interactions are naturally devel- 
oped and visualized in momentum representation, where 
the coupling of momentum scales can be interpreted via 
scattering theory [13]. However, softening from the RG 
is accompanied by induced non-locality, which obfuscates 
the nature of the softening in coordinate representation. 
Indeed their features in coordinate space have only been 
inferred rather than directly visualized. This has left 
open questions about the nature and consequences of the 
non- locality. For example, the elimination of short-range 
correlations with softening has been associated with the 
elimination of the repulsive core, but only by implication. 
Can it be directly visualized and how does this change 
SRC intuition? Is the long-range part unchanged? How 
is the RG flow toward universal low-momentum interac- 
tions manifested in coordinate representation? Is there 
a way to use the non-local low-momentum interactions 
with modern QMC techniques? 

To address these questions, in this paper we explore 
a simple "local projection" for two-body potentials that 
separates softened non-local interactions into a local ef- 
fective interaction plus a non-local residual interaction 
in each partial wave. This provides a natural visualiza- 
tion of the range-dependent features as the potential is 
evolved by RG methods. There is no unique projection 
at short distances, which gives us wide freedom to choose 
a form for the projection that reduces to local pion ex- 
change at large separations (r ^ 1 fm). In Section II, we 
present a particular projection and describe its applica- 
tion and properties. Projections are made of evolved sim- 
ilarity renormalization group (SRG) potentials [13] start- 
ing from Argonne vi^ [10] and from the Entem-Machleidt 



500-MeV-cutoff chiral EFT potential [11] are given in 
Section III. The resulting local potentials demonstrate 
the elimination of short-range repulsive cores and the 
flow toward universal low-momentum interactions. In 
Section IV, the SRG residual interaction is shown to be 
relatively weak and we test its inclusion in perturbation 
theory. We summarize and outline further implications 
and applications in Section V. 



II. BACKGROUND 

We consider functionals that take a two-body potential 
that is dependent on both coordinate indices and produce 
a local interaction multiplying a delta function: 



L[V{r,v')]=5ir-r')fiv). 



(1) 



We call L a local projection if it acts as an identity func- 
tional for already local potentials, 



L[L[y(r,r')]]=L[y(r,r')]. 



(2) 



Our original motivation was simply to visualize what is 
happening at different length scales in coordinate rep- 
resentation as the SRG (or another transformation) is 
applied to soften two-nucleon (and three-nucleon) inter- 
actions. However, such projections can also be used to 
decouple local and purely non-local components of the 
interaction for calculations. 

Perhaps the simplest possible choice for a local projec- 
tion is 

L[F(r, r')] = S{r - r') j dv" V(v, v") 

^6{v-v')V{v), (3) 

which sums at each r the interaction weight from the 
connected r" coordinates, essentially averaging over the 
non- locality. This should capture most of the effect of the 
potential on long- wavelength nucleons. It is clear that 
this functional acts as the identity for local potentials. 
That is, if V(r, r^ = V{y)8{y - r'), then 



V{v) = / dv" V{t)S{t - r") = V(r) . 



(4) 



Here and for the remaining discussion we drop the com- 
mon factor S{r — r') from the definition and work only 
with the function V(r) multiplying it. 

This instance of local projection is not useful for tensor 
forces or terms local in only the radial coordinate, such 
as spin-orbit terms of the form f{r)L • S. To deal with 
such terms, we extend the definition to apply in a coupled 
partial wave basis. A natural extension is: 

v\Z'{r) = jdn^pdv' l^*„(l^r)ll'™'(OrO^(r,r') 
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This leaves unchanged local operators such as those in 
the Argonne vi^ interaction, where 



Vl^.{ry) = Vl^.{r) 



im ,^Ar-r') 



(6) 



We define our partial wave expansion to have unit nor- 
malization, 

V{r,r') = Y,Y.^L{^)^i'm'{v')VlZ'{r,r') , (7) 

Im I'm' 

and to be symmetric between momentum and coordinate 
representations, 

F(k ,k') = E E ^™(k )Yi,m'ik')Vl;Z'ik, k') . (8) 

Im I'm.' 

Two-nucleon SRG interactions are usually computed 
directly in momentum representation [13] and it is there- 
fore convenient to express the projection Eq. (5) in that 
basis. This follows by inserting the Fourier transform of 
the momentum-space interaction and calculating the free 
coordinate integral analytically, yielding 



: Im , 



yooW= / fc^dfciK^O^oo (fc,0) 



(9) 



and (for ;,r >0) 



where 
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This projection provides a useful starting point for vi- 
sualizing what is happening at long and short distances 
as we make an RG evolution of the interaction. How- 
ever, it is not unique and we have no evidence that it 
is an optimal separation of a local part of the interac- 
tion for calculations. From this point forward, we will 
only consider potentials that are rotationally invariant, 
so that the tti, m' dependence is a Kronecker delta func- 
tion, which will be suppressed. 

To compare the relative strength of the local projection 
to the remaining non-local residual interaction (V — V), 
we will use the Fourier transform of the local projection, 

: im /^ J ,^ ^ I 2 T • /7 \ • /7/ \TT ^''^ 
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r^drMkr)ji,{k'r)V\Z,{r). (12) 



Expressing this momentum-space local projection di- 
rectly in terms of the original momentum-space non-local 
potential is helpful to understand how universal features 
of SRG-evolved interactions affect the projection and how 
non-local features affect expectation values. For even / 
and r, a direct expression is known: 
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q^dqI{l,l,l'q,k,k')flZ'{Q), (13) 



where 



^00^^0(9,0) 



if r = , 

fIZ'iQ) = \vZ'iQ,0) in = 0, (14) 

(Nw Jdq' q'-'VlZ'{q,q') otherwise, 

and /(/, /, V g, /c, k') is an integral of three spherical Bessel 
functions given in Ref. [14]. It is also numerically more 
stable to evaluate Eq. (13) instead of Eq. (12). For the 
figures presented here, we need: 
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where A = k'^ 



A:'^ B = e -k'^, and 



/\{q,k,k') = < 



1 \k-k'\<q<k^k' , 
^ \k — k^\ = q or q = k -\- k' 
otherwise. 



(15) 
(16) 

(17) 
(18) 



The expressions in Eqs. (13)-(18) provide insight by 
showing the different weightings of momentum regions 
for the S- and D- waves, as seen below. 



III. VISUALIZATION 

In this section, we use the local projection from Sec- 
tion II to construct simple visualizations of SRG-softened 
interactions [15, 16]. The SRG generates a continuous se- 
ries of unitary transformations labeled by a parameter 5, 

H{s) = U{s)H{s = 0)U'^{s) , (19) 

which can be implemented as a flow equation [16, 17] 



ds 



H{s) = [v{s),H{s)], 



(20) 
(21) 



In the present studies, we use only the most common 
choice of ri{s) for nuclear applications, 



V{s) = [Tre,,H{s) 



(22) 



which has demonstrated momentum decoupling proper- 
ties [13, 18-26]. We use the momentum decoupling scale 
A = 5^/^ (in units where h = c = m = 1 with nucleon 
mass m) to label the evolution. In typical applications, A 
ranges from oo (which is unevolved) to final values from 
1.5 to 2.0 fm"^ [13]. 

In Figs. 1 to 4, the local projections from Section II are 
applied to the original and SRG-evolved Argonne 'Ui8 [10] 



and Entem-Machleidt N^LO (500 MeV) [11] NN poten- 
tials in selected channels. These interactions are evolved 
in momentum representation from A = oo (initial) to 
A = 1.6 fm~ and the local projection is calculated from 
the momentum representation. In each series of plots, 
the initial projected potentials are maintained as dashed 
lines for comparison to the evolved version and the region 
between the evolved projections is shaded to highlight the 
approach (or non-approach) to universal form. 

The A = oo panel in each of these figures shows the 
original Argonne vis potential (up to small numerical 
artifacts from the use of finite momentum meshes) be- 
cause it is local and projects onto itself. The effects of 
the strong short-range core are evident in all channels, 
with the peaks in the S-waves being off scale at over 
2 GeV. The projection of the unevolved N^LO potential 
(which is not local) exhibits an overall oscillation with 
wavelength of about 2.5 fm, as would be expected from a 
momentum cutoff of 500 MeV. This cutoff is much lower 
than for Argonne ^'ig and therefore the potential is much 
softer, as shown by much weaker short-range contribu- 
tions, which peak in the S-waves at about 140 MeV for 
the singlet channel and 200 MeV for the triplet channel. 
In the coupled ^Si-^Di channel in Fig. lb, which in- 
cludes contributions from the tensor force, we again see 
very strong effects at short distances for Argonne vis and 
only relatively weak effects from the N^LO potential. 

In the S-wave channels for Argonne Vi8, Figs, la and 
2 show that the SRG significantly modifies short-range 
features first as A decreases, unlike in momentum rep- 
resentation where the initial modification occurs at all 
momenta. This is consistent with the modification tak- 
ing the form of a regulated delta function in coordinate 
space, which is primarily a uniform shift in momentum 
space below the cutoff scale. Such a shift has been pre- 
viously identified for the SRG [27, 28] and for the Mow/c 
RG [29]. We see a rapid and dramatic dissolution of 
the repulsive core and a relatively weak modification at 
longer ranges (r > 2 fm) once we have evolved below 
about 2 fm~ . That is, the core is largely removed before 
significant modification to the tail begins. 

This "core meltdown" is qualitatively consistent with 
previous intuition about SRG evolution, but to our 
knowledge this is the first quantification of the size of 
the modification at both short and long ranges. The 
N^LO evolution is less dramatic, but equally complete in 
the suppression of short-range repulsion. Both projected 
potentials are completely attractive by A = 2fm~^, ex- 
cept for small long-range oscillations. This is consistent 
with observed perturbative behavior in infinite nuclear 
matter, where the repulsion driving saturation originates 
with three-body forces [21, 30]. From this point on in the 
flow the potentials become increasingly alike at all dis- 
tance scales, reflecting the flow toward universal form 
previously observed in momentum representation [13]. 
Because the local projection samples low- momentum ma- 
trix elements with greater weight than those at high mo- 
mentum (see Eqs. (13)-(18)), low-momentum universal- 
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(b) Snapshots taken at selected A of the SRG evolution in the ^Si-'^Di channel 
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(c) Snapshots taken at selected A of the SRG evolution in the "^Di channel 



FIG. 1. (color online) SRG evolution of local projections for the deuteron NN channels, with the unevolved (A = oo) projected 
Argonne vig [10] and N^LO (500 MeV) [11] potentials shown as dashed lines. The region between the potentials is shaded to 
highlight the approach to a universal form in the "^Si and ^Si-^Di channels. The initial Argonne vig potential is local and 
therefore its unevolved local projection is unchanged by local projection. 
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FIG. 2. (color online) SRG evolution of the local projections for the ^So channel as in Fig. 1. 
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FIG. 3. (color online) SRG evolution of the local projections for the "^Po channel as in Fig. 1. 




FIG. 4. (color online) SRG evolution of the local projections for the ^D2 channel as in Fig. 1. 



ity translates to universality at all ranges in coordinate 
space for S- waves. 

In the ^Si-^Di coupled channel where the tensor force 
is active (see Fig. lb), we observe dramatic reduction of 
short-range strength as a result of the SRG evolution. 
The initial N^LO projection already is weak compared 
to Argonne v\^^ which is consistent with both the cen- 
tral suppression seen in the S-waves and the suppres- 
sion of the tensor contribution by the 500 MeV cutoff. A 
universal form at all distance scales is seen already by 
A = 2fm~ . In contrast, the projected potentials in the 
^Di channel in Fig. Ic are only equal for r > 2fm~ 
and this is true for all A. This contrasts with the S- 
wave universality, largely because of the stronger weight- 
ing of high-momentum matrix elements in the local pro- 
jection. The SRG evolution shows little modification of 
the N^LO potential at all A while there is significant but 
incomplete (compared to the N^LO flow) modification of 
the Argonne v\^ hard core. 

The evolution in the ^Pq channel (see Fig. 3) is qual- 
itatively the same as the S-waves: strong reduction ini- 
tially of short-range strength and later long-range modi- 
fications. However, it is less complete by the lowest A and 
the results are less universal at short range. The same 
is true but to a greater extent in the ^D2 channel (see 
Fig. 4). The reduced universality in the local projection 
is again a result of the stronger weighting of higher mo- 
menta when compared to the S-wave (see Eq. (13)). This 
results in the local projection matrix elements depending 
on diagonal matrix elements of the full interaction in the 
D-wave, peaking around k = k' = X, where universal- 



ity is incomplete. Note, however, that the effects of the 
residual short-range differences in matrix elements will 
be suppressed by the angular momentum barrier. 

By taking the Fourier transform of a local projection 
and subtracting it from the original momentum-space po- 
tential, we can examine the residual non-local compo- 
nent of the interaction. This is shown for the ^Sq chan- 
nel in Fig. 5 for Argonne vis and in Fig. 6 for N^LO . 
The A = oo plots for V — ]/ show the effect of the lo- 
cal projection on the initial potential: Argonne Vi8 is 
purely local and therefore V — V is identically zero, while 
the non-locality of N^LO is evident for k^k' > 2fm~^. 
The dominant source of this non-locality is the ultravio- 
let regulator used, which is of the form f{k)f{k') where 
f(k) = e~^^ /^ ^"^ with integer n and cutoff A. This 
is confirmed by observing the decrease in the non-local 
residual as A is increased. 

We note that N^LO is effectively local when at least 
one of k or k' is less than 1 fm~^. For both potentials, this 
region stays almost completely local during the full SRG 
evolution. Furthermore, the non-local component is rel- 
atively small elsewhere except near the large-momentum 
diagonal, which is decoupled from the low-energy physics. 
This suggests that the local projections can be more than 
a tool to visualize SRG evolved interactions, i.e., we may 
be able to treat SRG potentials as being local plus a 
small non-local part that can be handled in perturbation 
theory. 

In summary, the local projections show that the poten- 
tials are modified by the SRG from the inside out (i.e., 
short-range first) as the decoupling scale A is lowered. 
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FIG. 5. (color online) Contour plots showing SRG evolution of the interaction (V), its local projection(y), and the residual 
non-local interaction {V — V) multiplied by 5, starting from the Argonne i;i8 potential [10]. The initial potential (A = 00) is 
local, as verified in (a), while the SRG evolution generates increasing non- locality. 



The flow to universal form is cleanly seen in the S- waves, 
which become shallow and structureless attractive poten- 
tials. In the higher partial waves, a nearly universal form 
is reached except at the shortest ranges where the impact 
of the remaining discrepancies in matrix elements may be 
limited. Plots of the non-local residual interaction imply 
it is relatively weak in the low- momentum region, sug- 
gesting the applicability of perturbation theory, which 
we consider in the next section. 



IV. PERTURBATION THEORY 

In addition to providing a visualization of non-local NN 
forces, the local projection gives a well-deflned non-local 



residual interaction, V — V. This residual by construction 
has no further local piece that can be extracted with the 
local projection used to compute V: 



L[V -V]= L[V] -V = V-V = 0. 



(23) 



As we saw from Figs. 5 and 6, the low- momentum part 
of V{k^k') (i.e., for k^k' < 1 fm~ ) is almost completely 
local even after signiflcant SRG evolution (A < 2fm~ ). 
This suggests that much of the low-energy physics is con- 
tained in the local interaction. Now we can ask how quan- 
titatively we can reproduce low-energy observables using 
the local piece exactly with perturbative corrections from 
the non-local piece. 

We start by considering phase shifts, using the dis- 
torted wave born approximation (DWBA) [31] to expand 
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FIG. 6. (color online) Contour plots showing SRG evolution of the interaction (V), its local projection(y), and the residual 
non-local interaction {V — V) multiplied by 5, starting from the N^LO potential of Ref. [11]. Even though the initial interaction 
is generally non-local, it is a nearly local potential at very low momentum. 



about the local projection. Because we are working with 
the non-local residual interaction, the entire calculation 
is most convenient in momentum space. To proceed, we 
first compute the full momentum-space scattering wave- 
function for the local projection in a single partial wave 
(we suppress dependence on /): 



iMi^i) 



Ap 
pk 



S{p-k)^BpV- 
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Wp{k) 



, (24) 



where Bp = — tt""*^ tan(J^(p)) is real, Ap = {1 -\- mBp)~^ 
is a complex coefficient, and Wp{k) is a real regular func- 
tion. These ingredients are computed non-perturbatively 
using the method described in Ref. [32]. We apply the 
two-potential formula [31, 33] to expand the exact half- 
on-shell partial- wave T-matrix in terms of the local scat- 



tering function, keeping only the first-order correction: 

T(fc,fc') = {k\V\^+) + {^-\(V-V)\xt,) 

« {k\V\ijt,) + ii^kliV-VMt') . (25) 

where |x^) is the scattering wave function for the full 
potential and |?/^^) is the scattering function for the lo- 
cal projected potential. We show representative results 
in Figs. 7 through 10 comparing phase shifts in the ^Sq 
and ^D2 partial waves calculated from the full potential, 
the local projection, and the first-order DWBA-corrected 
local projection at several stages in an SRG evolution. 

Consider first the ^Sq phase shifts in Figs. 7 and 8. 
The Argonne '^is potential is local, so for A = oo the 
phase shifts from the projected potential are exact. In 
contrast, the phase shifts from the non-local unevolved 
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FIG. 7. (color online) Nucleon-nucleon phase shifts in the ^So channel at three stages of an SRG evolution (A = cxo, 3fm~^, and 
2fm~^) starting from Argonne vig as the initial potential. In each panel, the phase shifts computed from the on-shell T-matrix 
for the full interaction (solid) are compared to those from the local projection alone (dashed) and from the projection plus a 
first-order DWBA correction (dotted). The insets show an expanded view at low energy. 
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FIG. 8. (color online) Same as Fig. 7 except starting from Entem-Machleidt N^LO 500 MeV as the initial potential. 
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FIG. 9. (color online) Nucleon-nucleon phase shifts in the ^D2 channel at three stages of an SRG evolution (A = 00, 3fm~^, 
and 2fm~^) starting from Argonne vig as the initial potential. In each panel, the phase shifts computed from the on-shell 
T-matrix for the full interaction (solid) are compared to those from the local projection alone (dashed) and from the projection 
plus a first-order DWBA correction (dotted). The insets show an expanded view at low energy. 



N^LO potential are modified by the projection at all ener- 
gies. The first-order perturbative correction improves the 
agreement everywhere, but there remains significant dis- 
agreement above 100 MeV. As the Argonne vis potential 



is evolved, the local projection increasingly deviates from 
the exact phase shifts below 100 MeV, but the first-order 
correction restores agreement to the one percent level. 
We note that reconstructing the low-energy peak in this 
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FIG. 10. (color online) Same as Fig. 9 except starting from Entem-Machleidt N^LO 500 MeV as the initial potential. 



channel is non-trivial because of the fine tuning in the po- 
tential that generates the large S-wave scattering length. 
Above 100 MeV the phase shifts from the local potential 
are more attractive than the exact phase shifts, consis- 
tent with the removal of the hard core. The perturbative 
contribution increasingly overcorrects in this region. For 
the N^LO potential, the effect of evolving in A is much 
less pronounced and the final results for A = 2 fm~ and 
below are very similar to Argonne Vi8, as expected from 
the flow to a universal form in this channel. 

The story for the ^D2 channel in Figs. 9 and 10 is 
qualitatively the same, except that the agreement at the 
lowest energies (below 50 MeV) remains very good in all 
cases. The deviations at higher energies are more severe 
for the N^LO potential because of the non-local regula- 
tor. 
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FIG. 11. (color online) Deuteron binding energy as a function 
of SRG A, starting from the Argonne v\^ potential. Exact 
results [Ej) are compared to the local projection (Ej) and the 

first two orders of ordinary perturbation theory {E and 



In Figs. 11 and 12, we compare the deuteron binding 
energy from the full Argonne v\% and N^LO potentials, 
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FIG. 12. (color online) Same as Fig. 11 except starting from 
the N^LO 500 MeV potential. 



respectively, to the results from the local projection and 
the first two orders of perturbation theory. The devia- 
tions for the local projections are both significant, with 
similar results after evolution to A below 3 fm~^. The 
error reflects the fine-tuned cancellation between the ki- 
netic energy and the potential energy; this fine-tuning is 
not as completely preserved by the local projection. Per- 
turbation theory improves the energy at all A values, with 
the second-order correction at A = 1.6 fm~^ bringing the 
energy for both initial potentials within 100 keV. 

We can also look at the probability densities for the 
deuteron, which are shown in Figs. 13 and 14, and the 
deuteron momentum distributions (i.e., the momentum- 
space wave function squared), which are shown in Figs. 15 
and 16. While these are not observables, they can give 
insight into how the local projection modifies the nuclear 
wave functions. Short-range correlations for the initial 
potentials are associated with the "wound" at small r in 
coordinate space and the strong high-momentum com- 
ponents (e.g., for k > 2fm~ ) in momentum space. The 
SRG evolution fills in the wound and greatly suppresses 
the high-/c strength. The deviations in the wave func- 
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FIG. 13. (color online) The deuteron wave function squared 
calculated from the full interaction evolved to A = 2fm~^ 
starting from the Argonne viq potential compared to the re- 
sult using the local projection and first-order perturbation 
theory. 



FIG. 15. (color online) The deuteron momentum distribution 
calculated from the full interaction evolved to A = 2fm~^ 
starting from the Argonne vis potential compared to the re- 
sult using the local projection and first-order perturbation 
theory. 
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FIG. 14. (color online) Same as Fig. 15 except starting from 
the N^LO 500 MeV potential. 



FIG. 16. (color online) Same as Fig. 15 except starting from 
the N^LO 500 MeV potential. 



tions caused by the local projection are largely removed 
by first-order perturbation theory. The short-range cor- 
relations are also additionally suppressed by the local 
projection relative to the full evolved potential, but the 
qualitative interpretation is unchanged. 

Finally, we can use Weinberg eigenvalues [34] as a diag- 
nostic for changes in the NN potentials under local pro- 
jection [35, 36]. In Fig. 17, we show the largest attractive 
and repulsive eigenvalues at zero energy for the ^Sq chan- 
nel as a representative example. The value of the large 
attractive eigenvalue, which approaches but is always less 
than one, reflects the near bound state in this channel. 
This is a physical property so we do not expect dramatic 
changes with decreasing A. (Note: this eigenvalue would 
be equal to one for £^ = if there were a bound state 
precisely at zero energy.) The differences in eigenvalues 



from the local projections and full potentials are consis- 
tent with the variations observed above for the deuteron 
energy. The extreme decrease of the largest repulsive 
eigenvalue for Argonne vi^ quantifies the melting of the 
repulsive core as a result of SRG evolution. A signifi- 
cant but much smaller effect is seen for the N^LO poten- 
tial, which is already comparatively soft initially. In both 
cases, the local projections generate eigenvalues that are 
always smaller than those from the full potential. This 
observation persists at all other energies. Thus the local 
projections are always at least as perturbative as the full 
potentials. The vanishing of the eigenvalues just below 
A = 4 fm~ means that the negative of the potential fails 
to have a bound state [36]. This is consistent with the 
observed evolution of the potentials toward being purely 
attractive in this channel. 
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berg eigenvalues at energy ^ = in the ^So channel as a 
function of SRG A for Argonne vig and N^LO initial poten- 
tials. The local projection V is compared to the full evolved 
potential in each case. 



V. SUMMARY 

The renormalization group evolution of nucleon- 
nucleon potentials to decouple low momenta from high 
momenta leads to increasingly non-local interactions. 
The non-localities inhibit clean visualizations in coordi- 
nate space and therefore make it difficult to develop an 
intuitive picture of the softening at different ranges. To 
overcome this limitation, we have studied a local projec- 
tion (or non-local average) of the interaction. As we fol- 
low the flow of the projection, we can see directly that the 
repulsive core dissolves first, and only later in the evo- 
lution do longer-ranged features become modified. For 
the S- waves and some (but not all) other partial waves, 
the local projections from different initial potentials ap- 
proach a common form at all ranges (which is purely 
attractive in the S- waves). 

The observed pattern of non-locality agrees with con- 
ventional wisdom that the non-locality generated by RG 
evolution to lower resolution increases with the momen- 
tum transfer q [13]. For example, the dominant effect of 
evolution for off-diagonal matrix elements is the applica- 
tion of a non-local form factor: 

= e-(''+'^')'/(^'/'i')y,=oo(k,k'). (26) 
This contrasts with a Perey-Buck-type factorized non- 



local potential [37] such as 
y(r,r') = V^ocai[(r + r72)]- 



-ir 



'(7ri/2/3)3- ' ^^'> 

for which the non-locality /3 is independent of the mo- 
mentum transfer. The preservation of locality for low- 
momentum transfer is consistent with physical con- 
straints that the longest-ranged part of the interaction 
is governed by pion exchange. 

A local projection separates the full potential into a lo- 
cal piece and a purely non-local residual (this means that 
it is annihilated by a second projection). We find that 
low-energy observables calculated with the local projec- 
tion of evolved interactions alone are reasonably well re- 
produced, with greater deviations for fine-tuned observ- 
ables and at higher energies. Including corrections from 
the non-local residual potential in perturbation theory re- 
duces deviations of phase shifts to the few percent level 
up to about 100 MeV (after which there remain signifi- 
cant deviations) and can correct the deuteron binding to 
better than 100 keV. This implies that the bulk of the 
low-energy physics is actually determined by local inter- 
actions and the purely non-local corrections are pertur- 
bative. 

If the higher-energy contributions are sufficiently de- 
coupled, our results suggest a possible strategy for ap- 
plying low- momentum potentials in QMC calculations. 
Namely one can use existing technology with the local 
projection part of the interaction and then correct the re- 
sult in perturbation theory. While the accuracy achieved 
in the test cases examined here may not be competitive 
with direct use of local potentials such as Argonne Vi8, 
the particular choice of local projection considered has 
not been optimized in any way. Further tests includ- 
ing different choices for the projection are needed before 
judging the feasibility of this program. More importantly, 
one needs to include the contributions of many-body in- 
teractions, which is also of great current interest for visu- 
alizing the nature of induced three- and four-body forces 
effects during SRG evolution (including the effects of dif- 
ferent choices of SRG generators). Work on this is in 
progress using a hyperspherical representation. 
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